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Abstract 

We formulate quantum tunneling as a time-of-arrival problem: we de- 
termine the detection probability for particles passing through a barrier 
at a detector located a distance L from the tunneling region. For this pur- 
pose, we use a Positive-Operator- Valued-Measure (POVM) for the time- 
of-arrival determined in ,1 . This only depends on the initial state, the 
Hamiltonian and the location of the detector. The POVM above provides 
a well-defined probability density and an unambiguous interpretation of 
all quantities involved. We demonstrate that for a class of localized initial 
states, the detection probability allows for an identification of tunneling 
time with the classic phase time. We also establish limits to the definabil- 
ity of tunneling time. 

We then generalize these results to a sequential measurement set-up: 
the phase space properties of the particles are determined by an unsharp 
sampling before their attempt to cross the barrier. For such measurements 
the tunneling time is defined as a genuine observable. This allows us to 
construct a probability distribution for its values that is definable for all 
initial states and potentials. We also identify a regime, in which these 
probabilities correspond to a tunneling-time operator. 



1 Introduction 

This paper is a continuation of Ref ( 1 ), in which a procedure was sketched for 
the construction of a Positive-Operator- Valued Measure (POVM) for the time- 
of-arrival for a particle described by a Hamiltonian H . Here, we extend this 
POVM to cover the case of particles tunneling through a barrier. This procedure 
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allows us to provide an unambiguous determination for the tunneling-time as it 
can be measured in time-of-arrival type of measurements. 

Quantum tunneling refers to the escape of a particle from a region through a 
potential barrier, whose peak corresponds to an energy higher than that carried 
by the the particles. There are two important questions (relevant to experi- 
ments) that can be asked in this regard. The first is, how long does it take a 
particle to cross the barrier (i.e. what is the tunneling time?). The second is, 
what is the law that determines the rate of the particle's escape through the 
barrier? In this paper, we develop a formalism that provides an answer to these 
questions and we apply it to the first one. The issue of the decay probability 
will be taken up in [2]. 

The issue of tunneling time has received substantial attention in the litera- 
ture, especially after the 1980's-see the reviews [Ull]. The reason is that there is 
an abundance of candidates and a diversity of viewpoints with no clear consen- 
sus. There are roughly three classes of approaches: (i) Wave packet methods: 
one follows the particle's wave packet across the barrier and determines the 
tunneling time through a " delay in propagation" Ej , (ii) one defines suitable 
variables for the particle's paths and one obtains a probability distribution (or an 
average) for the transversal time corresponding to each path. These paths can 
be constructed either through path- integral methods [H El [5] , through Bohmian 
mechanics [3] , or through Wigner functions [TU] , and (iii) the use of an observ- 
able for time: this can take the form of an additional variable playing the role 
of a clock [TTl [12] , or of a formal time operator [13] . In general, these methods 
lead to inequivalent definitions and values for the tunneling time. 

1.1 Our approach 

The basic feature of our approach to this problem is its operational character. 
We identify the tunneling time by constructing probabilities for the outcome 
of specific measurements. We assume that the quantum system is prepared in 
an initial state '0(0), which is localized in a region on one side of a potential 
barrier that extends in a microscopic region. At the other side of the barrier 
and a macroscopic distance L away from it0, we place a particle detector, which 
records the arrival of particles. Using an external clock to keep track of the 
time t for the recorder's clicks, we construct a probability distribution p{t) for 
the time of arrival. The fact that the detector is a classical macroscopic object 
and that it lies at a macroscopic distance away from the barrier allows one 
to state (using classical language) that the detected particles must have passed 
through the barrier (quantum effects like a particle crossing the barrier and then 
backtracking are negligible). Hence, at the observational level, the probability 
p{t) contains all information about the temporal behavior for the ensemble of 
particles: all probabilistic quantities referring to tunneling can be reconstructed 
from it. 

With the considerations above, both problems of determining the tunneling 

^We explain in section 2.3 the sense in which we employ the word "macroscopic". 
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time and the escape probability as a function of time (see [5]) are mapped to the 
single problem of determining the time-of-arrival at the detector's location for 
an ensemble of particles described by the wave function i/jq at t = and evolving 
under a Hamiltonian with a potential term. To solve this problem, we elaborate 
on the result of [T], namely the construction of a Positive Operator Valued 
Measure (POVM) for the time-of-arrival for particles for a generic Hamiltonian 
H-see [14] and [15] for definition, properties and interpretation of POVMs. This 
POVM provides a unique determination of the probability distribution p(t) for 
the time-of-arrival. It is important to emphasize that by construction p(t) is 
linear with respect to the initial density matrix, positive-definite, normalized 
(when the alternative of non-detection is also taken into account) and a genuine 
density with respect to time. 

Since our results depend on the POVM for the time-of-arrival constructed in 
[1] , we review here the basic physical considerations involved in its construction. 
The technical aspect, namely the construction of this POVM for the problem 
at hand is undertaken in Sec. 2. 

The POVM of involves no structures other than the ones of standard 
quantum mechanics: the Hamiltonian, the initial state and the location of the 
recording device. It also involves a smearing function with respect to time, 
but we employ it in the regime in which the results are independent of such a 
choice. The first step in the derivation arises from the remark that the notion 
of arrival-time is well-defined when one considers histories for a physical system 
(both in classical and in quantum probability). We assume that the detector 
lies at X = L and that the initial state is localized in the region / = {x, x < L}. 
Moreover, we assume a discretization ^0,^1,^2, ■ ■ ■ ,tn of a time interval [0,r]. 
One asks at any instant U of time, whether the particle lies in region I or in 
region II = {x,x > L}. The set of all possible successive alternatives forms a 
Boolean algebra. The key point is that one can construct a subalgebra of events 
labeled by the time of first crossing (together with the event of no crossing), 
namely by the first instant ti that the particle is found in region //. This 
implies that propositions about the time-of-arrival have a well-defined algebraic 
structure, which is compatible with the Hilbert space description of quantum 
mechanics. The algebra of propositions for the time-of-arrival is a special case of 
the so-called spacetime coarse-grainings [TBI HZ] that have been studied within 
the consistent histories approach to quantum mechanics [HI [HI [2Ql [2Tj . 

The construction above takes place at the discrete-time level. One should 
then implement the continuum limit within the quantum mechanical formalism. 
The problem is that there is no proper continuous limit if one works at the level 
of probabilities (for the same mathematical reason that leads to the quantum 
Zeno effect). However, there is a proper continuum limit for this algebra if 
one works at the level of amplitudes. More specifically, one can implement the 
continuous limit at the level of the decoherence functional, an object introduced 
in the consistent histories approach. 

The decoherence functional is a hermitian, bilinear functional on the space 
of histories that contains all probability and phase information for the histories 



3 



of the systeno. The restriction of the decoherence functional to the algebra 
of propositions about the time-of-arrival effectively yields a hermitian function 
p(t, t') which is a density with respect to both of its arguments. 

The decoherence functional contains sufficient information for the construc- 
tion of POVMs for measurements that involve variables that refer to more than 
one instant of time. This has been established for sequential measurements 
[24] and for time-extended measurements [25]. In these cases one can compare 
the results to ones obtained from single-time quantum mechanics, but for the 
time-of-arrival, there is no analogous construction without the use of histories. 
Nonetheless, the method provides a definition of POVM for the time-of-arrival 
through a suitable smearing of the diagonal elements of the decoherence func- 
tional. For a free particle, this reproduces Kijowski 's POVM [26] in the appro- 
priate regime. 

The important point in the procedure above is that the POVM of [1] is 
valid for a generic Hamiltonian. The time parameter entering the POVM is 
the external Newtonian time and the identification of the time-of-arrival is done 
through purely kinematical arguments. Hence, this result can also be applied 
to the specific Hamiltonian operators that are relevant to tunneling. This is the 
content of Sec. 2. 

Summarizing, there are three basic features in our approach: a) the refor- 
mulation of tunneling as a time-of-arrival problem, b) the use of POVMs for the 
determination of the probabilities for the tunneling particles, and c) the basic 
ideas of the histories approach that enable us to construct a suitable POVM. 

1.2 Relation to other approaches 

There are some common points and some points of divergence with previous 
work on the tunneling time issue. Yamada has employed the decoherence func- 
tional showing that different definitions of tunneling time correspond to different 
definitions of the alternatives for the 'paths' considered in the definition [?7] . 
The construction of the decoherence functional is different from ours in one re- 
spect: the (coarse-grained) histories we consider refer to the paths' first crossing 
of the surface x = L, which lies a macroscopic distance away from the barrier. 
In [271 , the histories refer to the crossing of the harrier and the ambiguity in the 
definition of the tunneling time refiects the inability to decide which of all possi- 
ble spacetime coarse-grainings provides the true measure of tunneling time. This 
is due to the fact that quantum 'mechanical' paths may cross and then reenter 
the barrier region. In our case, this is not an issue. The detector is far away 
from the barrier region (at a macroscopic distance L) and the probability that 
a particle crossing x = L would ever backtrack to the barrier is practically zero. 
Another difference is that Yamada argues within the context of the decoherent 
histories programme that deals with closed systems [28]. While we employ the 
methods and (many) conceptual tools of consistent histories, our approach is 

^Alternatively, it can be viewed as a generating functional for all possible temporal corre- 
lation functions of the system [221 123) . 
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strictly operational within the Copenhagen interpretation. The probabilities 
we construct refer to measurement outcomes in a statistical ensemble. The de- 
coherence functional is only used as a mathematical object that allows us to 
construct a POVM and the particle crossing of the surface x = L is viewed as 
corresponding to an (irreversible) act of measurement by a device located there. 

The fact that the measurement of the particle takes place far away from the 
barrier region suggests that our results should be compatible with the asymp- 
totic analysis of wave packets. Indeed, as we shall see, our expression for tunnel- 
ing time (whenever this can be defined) corresponds to the classic Bohm-Wigner 
phase time [5]. However, the methodology is different: we do not identify time 
through the peak x{t) of the wave-packet (or through its center-of-mass), but 
the probability distribution for the detection time is obtained from a POVM 
that is defined for all possible initial states. Unlike time of detection, a sharp 
definition for the tunneling time is only possible for initial states characterized 
by a strong peak in their momentum distribution. However, the generality of our 
construction allows us to fully specify the limits in the definability of tunneling 
time. 

From the technical point of view, our approach has more in common with 
the second class of proposals we mentioned in the beginning: time being iden- 
tified at the kinematical level from the properties of 'paths'. In particular, the 
formalism bears substantial resemblance to the Feynman path integral deriva- 
tion of tunneling times of Sokolovski and Baskin T. However, our boundary 
conditions are different, and more importantly the probabilities we obtain arise 
from proper probability densities with respect to time. While the time-averaged 
quantities in [7] are linear with respect to a restricted propagator, such prop- 
agators appear in a quadratic form in our expression for the probability. It 
was argued extensively in [T] that this is necessary, in order to obtain a genuine 
probability density in a way that respects the convexity of the space of quantum 
states. The present construction also shares these properties and this implies 
that the issue of complex tunneling times does not arise. 

1.3 Our results 

In Sec. 3, we apply the POVM we constructed in Sec. 2 to a simple case of a 
particle in one-dimension. We consider a potential barrier V{x), which takes 
non-zero values only in a bounded (microscopic) region of width d around a; = 0. 
A wave packet approaches the barrier from the negative real axis, while the 
detection of the particle takes place at x = L >> a. 

We assume that the initial wave -function is well localized in position and 
in momentum (e.g. a coherent state). In addition, we require that a/ko << 1, 
where kg the mean momentum and a the momentum spread of the initial state. 
We then find (for a rather general regime for the values of the parameters char- 
acterizing the system) that the probability distribution for the time of arrival is 
sharply peaked around a value t,n- 

From this probability distribution, we identify the delay due to the presence 
of the barrier as the difference between the time tm and the time it would 
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take a classical particle of momentum kg to travel from the center of the initial 
wavepacket to the location of the detector. This 'delay time' equals 



where Tk is the transmission amplitude corresponding to the potential V{x). 
The delay time may be negative; the tunneling time is obtained as the sum of 
td with the time it would take the particle to cross the barrier: it coincides with 
the classic phase time. 

Note that the delay and tunneling times defined this way are not observables 
of the system: they cannot be defined for a generic initial state, but only for 
states well localized in momentum and for values of the parameters that lead to 
a probability distribution p{t) characterized by a sharp peak. In this case, one 
can use classical arguments for their definition. However, if either the initial 
distribution has a substantial momentum spread, or if p{t) exhibits a more 
complex structure, there is no unambiguous way to define tunneling time. While 
the value for the time-of-arrival is a genuine observable (a random variable on 
the sample space of the POVM), the tunneling time as we define it here requires 
the knowledge of the corresponding timc-of-arrival for a free particle: and this 
cannot be defined, unless the initial value of momentum is knowrd. 

Hence, in this approach the tunneling time is a parameter of the detection 
probability. It can only be identified for specific initial states, and not for any 
state, because its definition involves a correspondence argument to classical 
physics. 

In Sec. 5, we propose a generalization of the results above that leads to 
a definition of the tunneling time as a genuine random variable. The idea is 
to consider a sequential measurement set-up: the phase space properties of 
the particle are determined through an unsharp phase space sampling before 
this attempts to cross the barrier, and then the time-of-arrival for the particles 
that crossed the barrier is measured. The sample space corresponding to such 
sequential measurements accommodates the definition of the tunneling time as 
a genuine quantum observable and it allows us to construct a marginal POVM 
that provides its probabilities for a generic initial state. In a specific regime, this 
POVM becomes independent of the details of the first measurement: as such it 
defines an ideal probability distribution for the delay and the tunneling times: 
this distribution suggests a definition for a delay-time and for a tunneling-time 
operator. 

•^Strictly speaking, the above definition of tunneling time involves counterfactual reasoning. 
However, in the operational setting we consider here this is not a problem, as long as we keep 
in mind that the tunneling time (whenever it can be defined) is a 'property' of the ensemble 
of detected particles and not of any individual one. 




(1. 1) 
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2 The general probability measure 



In this section, we review the construction of the POVM for the time-of- arrival 
in [1], and we extend it for the case relevant to tunneling. 

2.1 The histories formahsm 

The POVM of [T] is constructed using some notions of quantum mechanical 
histories, as they appear in the consistent histories approach to quantum theory 
of Griffiths TB], Omnes (T^, Gell-Mann and Hartle Hi HI]. We should note 
however that these objects are used in the present context differently, namely 
in an operational approach to quantum theory-see [23[ 124] 

A history intuitively corresponds to a sequence of properties of the physical 
system at successive instants of time. A discrete-time history a is then repre- 
sented by a string Ai , Pt2 ; ■ • ■ Pt„ of projectors, each labeled by an instant of 
time. From them, one can construct the class operator 

Co, = UHh)Pt,U{h)...U^{tn)Pt^Uitn) (2. 1) 

where U{s) — e^"^^" is the time-evolution operator. A candidate probability for 
the realisation of this history is 

p{a) - Tr (cIp^C^) , (2. 2) 

where po is the density matrix describing the system at time i = 0. 

However, the expression above does not define a probability measure in the 
space of all histories, because the Kolmogorov additivity condition cannot be 
satisfied: if a and (3 are exclusive histories, and aM (3 denotes their conjunction 
as propositions, then it is not true that 

p{ay (3)^p{a)+p{P). (2.3) 

The histories formulation of quantum mechanics does not, therefore, enjoy the 
status of a genuine probability theory on the space of all histories. 

However, an additive probability measure is definable, when we restrict to 
particular sets of histories. These are called consistent sets. They are more 
conveniently defined through the introduction of a new object: the decoherence 
functional. This is a complex-valued function of a pair of histories given by 

dia, P) = Tr (clpoCa^y (2.4) 

A set of exclusive and exhaustive alternatives is called consistent, if for all pairs 
of different histories a and /3, we have Re c?(a, (3) — 0. In this case, one can use 
equation (|2. 2p to assign a probability measure to this set. 
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2.2 Time-of-arrival histories 



Using the histories formahsm we construct a decoherence functional for time-of- 
arrival histories with N time steps ti,t2,- ■ . ijv (discrete-time). The reason for 
this construction is that the decoherence functional has a good continuous time 
limit (unlike the probabilities for histories). 

We consider a particle in one dimension for concreteness, even though the 
results obtained here only use abstract Hilbert space operators and hold more 
generally. We split the line into the interval {—oo,L] and the interval [L,oo). 
Let P_ and P+ be the corresponding projectors. Our aim is to identify histories 
that correspond to the statement that the particle crossed from the — region 
to the -I- region during a particular time step. If we assume that at t = the 
particle lies at the — region then it is easy to verify that the history 

Q!m := (-P- , ii; -P-7 ^2; • • • , -P-, -P+, im+i; 1, • • • 1, iw) (2. 5) 

corresponds to the proposition that the particle crossed x — L for the 
first time between the m-th and the m -I- 1-th time step. The sequence a = 
(F_, ii; F_, . . . , P^.tm] ■ ■ ■ /"_,%) corresponds to the proposition that the 
particle did not cross x — L within the n- time steps. 

The set of histories am together with a is exhaustive and exclusive (a sub- 
lattice of the lattice of history propositions)-see also [161 [29] . The decoherence 
functional is then defined on this set of histories: it is a hermitian bilinear func- 
tional on a sample space consisting of the points (ti, . . . ,t„) together with the 
point N corresponding to no crossing 

d{tn,tm) = d(Q;„,a,„) (2. 6) 

d{N,tn) = d(a,a„) (2. 7) 

d{N,N) = d{a,a). (2. 8) 

We next consider two discretisations {to = 0, ii, i2, ■ • ■ ^tv — and {<g = 
0,t[,t'2, ...fj^, = T} of the time interval [0,T] with time-step St = T/N, and 
5t' = T/N'. We construct the decoherence functional d{[t,t + 5t], [t',t' + St']), 
where n = tN/T and m = t'N'/T. This reads 

di[t, t + St], [t', t' + St']) = Tr (po[e'^''*'P-]"e*^**'P+ 

We then take the limit N, N' ^ oo, while keeping t and t' fixed. Assuming that 
po lies within the range of P_, i.e. P_/5o-P- = po we obtain 

d{[t,t + St],[t',t' + St']) ^ StSt'Tr {e'"^^' P+Hp^CtPQCl,P-HP+^ ,{2. 10) 

where Ct = lim„^oo(i^-e-*'^*/"P_)". Writing 

p{t,t') = Tr [e'"^^' ^''^ P+Hp^ClpQCtP-HP+^ (2. 11) 



8 



we see that the decoherence functional defines a complex-valued density on 
[0,T] X [0,T]. The additivity of the decoherence functional (which reflects the 
additivity of quantum mechanical amplitudes) allows us to obtain a continuum 
limit, something that could not be done if we worked at the level of probabilities. 

2.3 The tunneling Hamiltonian 

For the simple case of a particle at a line with Hamiltonian H = + V{x), 
where the potential V{x) is bounded from below, we employ a result in [l6l[30] 
that the restricted propagator Ct is obtained from the Hamiltonian H by im- 
posing Dirichlet boundary conditions at a: = L. If we also denote by GQ{x,x'\t) 
the full propagator in the position basis (corresponding to e~*^*), we obtain 

^(^'*') = -^d^{Ct.^oY{L)d^{CtiJa){L)G^{L,L\t-t') (2.12) 

where po = |'0o)('/'o|, with -00 having support for x < L. 

We now specialize to a case relevant for tunneling. We assume that the 
potential is short-range: it is significantly different from zero only in a region of 
width d around a; = 0. The distance L is macroscopic, while a is microscopic. 
This means that in the neighborhood of a; = L the propagator is effectively that 
of a free particle. Hence, we can substitute Go{L, L\t' — t) in Eq. (|2. 12p with 
the corresponding expression for the free particle 

G(L,L|i)=(— ) . (2.13) 

The considerations above also specify the range of values for L that are 
relevant to our problem. The first condition on L is that the propagator may 
be substituted by that of the free particle, as in Eq. (|2. ISp . The second is 
that L is sufficiently far away from the tunneling region so that the probability 
amplitude of a particle backtracking to the barrier region from L is practically 
zero. Physically one expects that this is the case for all initial states -00 for 
which the position spread l^q{t) remains at all times much smaller than L. 
Clearly, with the considerations above it is not necessary that L is a macroscopic 
distance in the literal sense of the word: the requirement that L be macroscopic 
is a sufficient but not a necessary condition. 

We next consider the Hamiltonian Hd that is obtained from the original 
Hamiltonian H by imposing Dirichlet boundary conditions at x = L. We dis- 
tinguish two cases: (i) if x takes value in the half-line, the spectrum of Hd 
is expected to be discrete; (ii) if x takes values in the full real axis, at least 
the positive energy spectrum will be continuous. (We restrict to Hamiltonians 
having this property.) Either way, for a; >> a, V{x) = and the solution of the 
Schrodinger equation HdiPe{x) — EiPe{x) with Dirichlet boundary conditions 
is proportional to sinfc(L — a), where k — {2MEy/^. We choose to label the 
eigenstates of Hd by fc, namely we write \k)u as a solution to the equation 

H\k)n^^\k)D, (2. 14) 
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with Dirichlet boundary conditions. 

Normalizing \k)£, so that 

D{k\k')D^5{k,k'), (2.15) 

(and similarly in the discrete-spectrum case) we write 

{x\k)D ^ Dusi\ik{L- x), (2.16) 

where the form of the normalization factor is specified the Hamiltonian's 
(generalized) eigenstates. 

For the study of tunneling, we assume that the initial state of the system 
has support only in the positive energy spectrum of H . Hence, 

(xlCtl^o) =^e-^'='*/2*^-DfcSinfc(L-a:)cfe, (2. 17) 

k 

where Ck = _D(fc|'0o) a-nd denotes the integration with respect to the spectral 
measure of Hd- Substituting into Eq. (|2. 12p . we obtain 

p(t,t')^ , ^ DkDtckctkk'e-' "^'^^'^*' . (2. 18) 



2.4 Construction of the POVM 

The decoherence functional contains sufficient information for the construction 
of POVMs for the probabilities of measurement outcomes for magnitudes that 
have an explicit time-dependence. In particular, the probabilities for the mea- 
surement outcomes for single-time, sequential and extended-in-time measure- 
ments (obtained through the standard formalism) can be identified with suitable 
diagonal elements of the decoherence functional-see [24l[25]. In other words, 
one can define POVMs by suitable smearing of the decoherence functional and 
in the cases above, these POVMs coincide with ones obtained from the stan- 
dard methods in quantum measurement theory. In the case of the time-of-arrival 
there is no analogous expression obtained from standard methods. However, the 
smeared form of the decoherence functional still defines a POVM, and the main 
assumption in this paper is that this POVM yields the correct probabilities. 

With this assumption, we obtain the following probability density for the 
time-interval [0, T] 



p^{t)^ ds ds'^F{t;s)^F{t;7)p{s,s'), (2.19) 



here /t(s, s') is a family of smeared delta functions /t(s, s') characterized by 
the parameter r. The functions /r satisfy the following property 

T 

dsr{s,s')^X%T]{s'), (2.20) 
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where X[o,t] is the characteristic function of the interval [0, T]: X[o,t]{s) = 1 if 
s G [0,T], and X[o,t](s) = otherwise. The functions /"^ incorporate specific 
features of the instrument that records particles crossing the surface x = L. 

Essentially, the key assumption in our approach (stated above) is that the 
functions /"^ appearing in the definition of (|2. 19p are analogous to the smearing 
functions that appear in the definition of POVMs for usual observables (i.e. ones 
other than the time of arrival). In [1], we showed that this assumption leads for 
the case of free particles to Kijowski's POVM 

The decoherence functional satisfies an hermiticity condition p(s, s') = p* (s', s) , 
which together with the positivity condition for its diagonal elements 

ds I p{s,s') > (2. 21) 



guarantees that p'^{t) is positive-definite for all values of t. 
The density ()2. 19p is linear with respect to the initial density matrix. To- 
gether with the probability of no-detection 

p^{N) = 1- [ dsp^'is) (2. 22) 

they define a POVM fl on the space [0,T] U {N}. This POVM describes the 
time of detection of a particle by an instrument located aX x — L. 

In this paper, we will be interested in taking T — > oo, i.e. taking t G [0, oo). 
It is convenient to work with Gaussian smearing functions 

r(s,s') = ^e-^. (2.23) 



However, these Gaussians are smeared delta-functions with respect to the whole 
real axis and not with respect to the time-interval [0, oo). To remedy this prob- 
lem, we note that by Eq. (|2. ISp . p(— s, — s') = p*(s, s') — p{s' , s) and that the 
probability (|2. 19p is symmetric to an exchange of s and s' . We also note that 
the contribution of terms that mix positive and negative s are significant only 
at times |t| of order t. Hence, the probability (12. 19p with the Gaussian ()2. 23p 
is substituted in place of /"^ and the integration limits taken from — oo to cxa, is 
twice the probability density that is defined with an integral over the positive 
half-axis. Hence, the use of Gaussian smearing functions only involves dividing 
p'^{t) in (|2. 19p by a factor of 2 (for times t» t). 

Inserting (|2. 231 into (|2. 19p . we change variables to m = ^(s -I- s') and v = 
s — s' noting that 

s) Vr s') = fr{u - t)e-f^. (2.24) 

We substitute in the integration fr{u — t) with a delta function 6(u — t). We 
then obtain 
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where 



-V /»r —lev 



dv- 



2^/7 



^^ e'y'/^[cos{2eTy)+sin{2eTy)]^ ^ 26) 



Vy 



At the hmit of er >> 1, i.e. if the detection time is much larger than e ^ 
e 2 [cos(2er2/) + sin(2eTj/)] 



dy 



Vy 



Hence, i?(e) = 2-y/7r/e. It follows that 
1 

M 

kk' 



Pit) = 



fcfc' 



Vfc2 + fc'2 



' — 333- 



(2. 27) 



(2. 28) 



The probability for the time-of-arrival then becomes independent of the pa- 
rameter r, and it is expressed solely in terms of the system's Hamiltonian, the 
initial state and the value of L. 

Eq. (|2. 28p is simplified if the spread Afc of the initial state |'0o) ( k = 
V 2MHd) is much smaller than the corresponding mean value fc: in this case, 
fc2 + fc'2 ~ 2fcfc', hence 



Pit) 



kCk\ e 



-ik^t/2M 



(2. 29) 



It was shown in T] that for the test case of a free particle (, in which Dj. 
(27r)^^/^) the probability distribution above reproduces the one of Kijowski [2| 



3 The detection probability 

In this section, we use the probability density (|2. 28p in a specific context that 
allows us to determine a magnitude that corresponds to the time the particle 
spends in the forbidden region. In effect, we identify tunneling-time by the 
delay caused by the presence of the barrier to the particles' time-of-arrival (see 
Sec. 4). This turns out to be the same definition as the one employed in the 
methods involving the wave packet analysis. However, we do not identify any 
specific features of the wave-packet (these objects have no natural probabilistic 
or operational interpretation in quantum mechanics), but we work directly at 
the level of measurement outcomes, namely the probability distribution for the 
time-of-arrival. 

We consider the simplest possible example of a particle tunneling through a 
potential barrier. We assume that the potential V{x) > takes non-zero values 
in a region of width d around x = 0. Let Vq be the maximum value of this 

*This condition is valid if the mean energy of the initial state is much larger than the 
energy uncertainty, and it is accurate for all times t >> t. 
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potential. In classical mechanics no particle with energy E < Vq can cross the 
barrier, hence the probability of detection at x = L is zero at all times. We next 
consider this problem in quantum theory. 

Eq. (|2. 28p involves the eigenstates of the Hamiltonian with Dirichlet bound- 
ary conditions a,t x = L. Since x € (— oo,cxd), the spectrum of the Dirichlet 
Hamiltonian is continuous. The summation over k is then substituted by an 
integral dk. 

The first step is to construct the generalized eigenstates of the Hamiltonian 
with Dirichlet boundary conditions. To do so, we first study the solutions to 
the Schrodinger equation 

- ^dluix) + VixMx) = ^u{x). (3. 1) 

There are two linearly independent solutions for each value of fc. It will nec- 
essary to construct an orthonormal basis of generalized eigenstates from these 
solutions. We pick one class of solutions u^(x) that correspond to a particle 
propagating from — oo and scattering on the potential 

u+(x) / ^'^ ("^"^ + Rte~^'^) X < -d/2 

where and is the reflection and transmission coefficient respectively, 
while is a normalization factor so that J dxu^{x)u^{x) = d{k — k'). Let 
be a normalized linearly independent solution that satisfies / dxu^ {x)u'j^ (x) = 
0. Its form will be the following 

u-(x) I (^^^'"'"^ + ^^"^'') ^ < -'^^^ (3 3) 

UfcW-| A- {e-''''= + R~e''''=) x > d/2 

Note that there is no reason for w^T to have a physical interpretation in terms of 
left-moving particles, and the labels T^T, are chosen for convenience: they do 
not correspond to a transmission and reflection coefficient of any short. We also 
note that the coefficients in u'^,u'i^ are not independent. For any two solutions 
ip, (f> to the Schrodinger equation with the same energy, the Wronskian ij/ (f) ~ (f)' ip 
must be x-independent. This yields the following conditions 

= - SkRt (3. 4) 

Sk^f+Rl+T^Rt (3.5) 

\Ti? + \Rt?^^ (3.6) 

K? + \Rl? = l + \Sk?- (3. 7) 

To impose the Dirichlet boundary conditions on these solutions, we take a 
linear combination Vk{x) of u'^ix) and uj^ix) and require that Vk{L) = 0. This 
yields 

vk{x) = Cu K-(l + R^e'^'^'^)u+{x) - A+T+e'^'^'^u^ix)] , (3. 8) 
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where 

Ck = , ^ — (3. 9) 

is a normalization constant chosen so that J dxvk{x)vk' (x)dx — 6{k — k'). 
For X > d/2, we obtain 

Vkix) = -2iCkA^A+T+e'''^ sin k{L " x). (3. 10) 

Hence, 

Dk = -2tCkA^A+T+e^''^ (3. 11) 

We now consider a Gaussian initial state tpo centered around xq < -'d/2 and 
having mean momentum ko > 

M^) = (^^-^"'^"^ (3- 12) 

where S is the spread in position and we assume that 6 << |a::o + d/2| so that the 
initial state does not overlap with the region where the potential is non-zero. 
In this region, 

Vk{x) = CkA-A+ [(1 + iJ-e^*'^^ - T+5fee2''=^)e^'=^ 

+ {R+R^e^'''^ - T+T-e2^'=^)e-*'=^] . (3. 13) 

The coefficients = Dikli^a) £^re then given by 
Ck = C\A-At [1 + (R- T+S,)e--'^'^^] ^--L_e-^-"«(^-+-')(3. 14) 

where we set a = (26)^^ the momentum spread. The assumption that tr/fco << 
1 allowed us to drop a term of order e"*^"/^'^ . 

The probability for the time-of-arrival at x = i is then given by p{t) = 
|z(i)P, where 

^it) = f dk B, (^^e-^ — ('=-»)/Ae--^*/-. (3. 15) 

In (13. 15p we defined 

Bu = -2^V2^\C,\M-,?\At? [l + [Rk - nSu)e-^^''^] T+. (3. 16) 

Since a/feo << 1, we can expand i?^, around its value at fc = /cq- As a term y/k 
also appears in the integral outside the exponential, we expand together 



14 



where 



2ko 



A. 



dlog\Bk 
dk 

darg[Bk] 
dk 



k—kn 



k—k() 



(3. 18) 
(3. 19) 



Within the same approximation, we take the limits of integration in Eq. p. ISp 
from —oo to oo. We then obtain 



1 



dke' 



UI (27rcr2)i/4 



(3. 20) 



The expression above involves a standard Gaussian integral. Its evaluation gives 

Z(t) - B, c-'klt/2M+rkoL [k^ (8^^)1 /4 

'° V4Af,/TT2lto 



X exp 



M. 

M 



l + 2ita^/M 



(3. 21) 



Hence, 



87rcr2 



X exp ■ 



2fcgaVM^ 
H-4t2cr4/M2 



4Af y l + 4i2cr4/Af2 



(l + ^&o'^ /'^ojt ; 



(3. 22) 



This expression is the probability distribution for the time-of-arrival, as it 
would be measured by a device located at distance L from the barrier. In the 
following section, we analyze its properties: in particular, we identify the delay 
caused by the presence of the barrier. 

4 Delay-time and tunneling time 
4.1 The identification of delay time 

For a sufficiently monochromatic wave packet (cr/fco — ^ 0), we assume that 
CfeoO'^/^o << 1, hence Eq. (|3. 22p yields 



Pit) = |z(t)p 



87rCT2 



/ 

4Af Y 1 + 4iVVM^ 



X exp ■ 



1 + 4t2cr4/Af2 



M{\xo\+L + Xk„) 



(4. 1) 
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The term 1 + At'^a'^/M^ corresponds to the spread in the particle's wave 
function due to time evolution. Since we want a configuration in which the 
determination of time is as sharp as possible, we assume that the value of a is 
so small that this spread is negligible at the time tm — ^^^^^"^^^^"'"'''''"^ , namely 
that t'^cr'^/M << 1. Then we obtain 



Pit) 



fcp 
4M 



TTCT^ exp ■ 



AP 



t- 



M(|xo| + i + AfeJ 



fco 



.(.4. 2) 



Then t — tm is a. sharp peak for the mean value of the time-of-detection. A 
classical particle (or in quantum theory a narrow wavepacket) that starts from 
Xq with momentum fco in absence of the potential barrier will arrive at x = L at 
(average) time t^ = M-l^^j^. Hence, the barrier causes a 'delay' td = t^ — to 
to the time-of-arrival (of the particles that are not reflected) 



td = MXkJko. (4. 3) 

The presence of the barrier has increased the effective length that has to be 
traversed by the particle by a factor of Afe^. In fact, Afeg may be negative: 
the time it takes the particle to cross the forbidden region of the barrier is 
ttun = M{Xko + dko)/ko, where dk^ = X2[ko) - xi(fco) > 0, where a;i,2(fco) are 
the points that determine the forbidden region: they are respectively the lowest- 
and highest-valued solutions of the equation — V{x). The total tunneling 
time has to be positive, but it is not necessary that it is larger than the time 
Mdkg/ko that the forbidden region is traversed by a classical free particle. 

We next calculate Xko in terms of the absorbtion and reflection coefficients 
corresponding to the potential V{x). From Eq. (|3. 16[) we see that the only 
term contributing to a phase in Bk is the product [l + fke~^^''^] , where 
/fc = {R- - f+Sk). We then obtain 

aiogr+\ , ^ /L-2iL/fe„ 



A^o = ( — ^ ) + ^-f^^— (4. 4) 



9k J k ,-„ 1 + /fee 



The second term in the right-hand-side of (|4. 4p oscillates very fast with L, 
because L is much larger than the de-Broglie wavelength 27r/fco of the particle. 
These oscillations are an artifact of our modeling the detection process by a 
crossing of the sharply defined surface x — L. In a realistic detection scheme 
the particle detection cannot take place with an accuracy grater than their 
de Broglie wavelength. For this reason, we can formally average L within a 
region of size I << L. Indeed, using a Gaussian smearing function p{L) = 
{'Kp)~^/^e~^-'"~^°'> we obtain a suppression factor of order e^'^o' << 1 for 
the oscillating terms. 

Hence, the effective tunneling time is 

fco fco V dk J^^, 



16 



i.e. we recover the expression for the Bohm-Wigner phase time [S]. It is im- 
portant to emphasize that this derivation did not employ any characteristics of 
the wave-packets (e.g. the trajectory followed by their peak, or their 'center- 
of-mass'). It is a natural operational definition at the level of the probability 
density that corresponds to the measurement outcomes. 

Note that a precise treatment involves smearing the probability function p{t) 
of (|3. 22[) . The only L-dependent objects that appear in this equation are the 
term Bkg and the Gaussian exponential. If i >> ^, the effect of smearing is to 
substitute L by the mean value Lq: the expression is not affected. The effect of 
smearing on Bkg is to suppress the oscillations; it leads to an effective expression 

Note that to a first (very rough) approximation, l^^^^l can be taken equal 
to (27r)~^/^, i.e. the value taken if the contribution of the region with no zero 
potential is considered to be negligible. Then 

(4- 7) 



V27r 

Before continuing, wc summarize the approximations involved in the results 
we obtained in this section. Eq. ()3. 15|) only involves the assumption that 
a/ko << 1. Eq. (j3. 22p involves the additional assumption that the function 
log Bk varies slowly around k — ko so that it is sufficient to keep the first order 

in its Taylor expansion. This approximation amounts to the condition |-q7^ — 

''0 

■g^jcr << 1. Eq. (|4. ip involves the additional assumption that ^fcnCr^/fco << 1- 
Finally, Eq. (14. 2|1 involves the assumption that t^^^a'^/M << 1. This implies 
that L cannot be too large, because the spread of the wave function due to 
the free propagation will induce a large uncertainty in the determination of 
tunneling time. 

4.2 Special cases 

Parity-invariant potentials. The expression for the mode functions and 
for Bk simplifies greatly if the potential is invariant under parity, namely if 
V{x) = V{—x). This implies that the eigenstate w^(x) can be identified with 
the parity transform of u^{x). Hence Sk = 0, = r^7 ■— Tk, = := Rk 
and A'^ = := Ak- We then obtain. 



\Ak\ 

|l + i?fee2^'=^|2 + |Tfc|^ 



= -2»^^ M , o LkL . + Rke-'^'']n (4. 8) 
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The square potential barrier Wc apply our results to the simplest example 
of a square potential barrier: V{x) = Vq for x G [—d/2,d/2]. Defining 7/5 = 
\/2MVo — k'^, we obtain the following values for the coefficients Tk, Rk 

J, ^ 2fc^_^fc^ 2fc7fc[2fc7fc cosh7fc(i - i(7^ - fc^) sinh7fcd] 

7/ 4^272 + (72 + fc2)sinh2 7fed ^ ■ ^ 

R ic-ikd + fc2)[2fc7Cosh7fcd - 2(72 - fc2) sinh7fcd] 

4/e272 + (72 + /e2)sinh2 7fcd ' ^ ' ^ 

There arc two limits, in which the results are particularly simple. The limit 
of a long barrier jkd >> 1, for which 



Tk ^ e-^^'e--^"" , f^^':... [2fc7fc - i(7^ - k')] (4.11) 



hi+k^r' 

-ikd -i.il -k'') + ik^k 



In this limit, the parameter Xkg is 

Xko = -d+—, (4.13) 

7feo 

i.e. it takes negative values (since jkod >> 1). The tunneling time is therefore 

. _ 2M 

The other limit is that of the delta function (very short) barrier. It is ob- 
tained by letting Vq ^ 00 and d such that Vod is a constant (we denote 
this constant as n/M). At this limit, 7^6? ~ VkcI and 

Hence, 

Afeo = TT^- (4. 16) 



Since d = Q the tunneling time is ttun = AfiSTrT'- 

"'"'^ 



4.3 Comments 

4.3.1 Domain of validity 

It is important to emphasize that our identification of a tunneling time ttun 
relies on the fact that the probability of detection has a unique sharp maximum 
at a specific moment of time. This is only possible for specific initial states. For 
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example, it is easy to demonstrate that a superposition of Gaussians centered 
at different values of momentum will lead to a probability distribution with an 
oscillating behavior. While there is still a mean detection time, we cannot read 
from it a time delay for the particle, because the momentum uncertainty does 
not allow one to specify uniquely a corresponding time for free particle evolution. 
Hence, the tunneling time is not a proper observable (i.e. a random variable 
on the sample space upon which the POVM is defined) in our description: it is 
only a parameter that appears in the detection probability for a class of initial 
states, which has an intuitive interpretation in terms of classical concepts. 

The fact that the concept of tunneling time has a restricted domain of valid- 
ity is highlighted by another point. We saw that for a long square potential the 

tunneling time equals ttun — . If d is very large, the condition 7fc„d >> 1 

"yfcQ "'0 

can be satisfied even if jka takes very small values, i.e. if the particle's mean 

energy ^ is very close to Vq. Hence, it is in principle possible to construct 
configurations, in which ttun is arbitrarily small: the effective 'velocity' d/ttun 
in the crossing of the barrier is then super-luminal. This is a well known effect 
in tunneling (the Hartman effect fSP). A full treatment in the present context 
involves the consideration of relativistic systems-this we will undertake in future 
work. Here, we only note that the regime of very large values for d/ttun, (very 
small values for jko ) is one for which the approximation involved in Eq. (|3. 17p 
fails. The tunneling probability increases rapidly in this regime and one would 
have to include further terms in the expansion of log Bk, which would lead to 
a substantially deformed probability distribution p(t) with no clear peak. The 
definition of ttun would then be highly problematic, and so would be the notion 
of a mean velocity in the tunneling region. 

4.3.2 Uncertainty in the specification of tunneling time 

The uncertainty in the determination of the peak in the probability distribution 
(j4. 2|) is In order for the delay time to be distinguishable (if we 

ignore all other sources of uncertainty) it is necessary that cr jAfcol >> 1. In 
order for the tunneling time to be distinguishable, it is also necessary to take 
into account the uncertainty in the quantity ^J"" . To leading order in a this 
equals akgCr, where 



The overall uncertainty in the determination of the tunneling time tt 
M{\ko + dko )/ko is of the order 



This expression is bounded from below by Mk^aka- Hence, a necessary 
condition for tunneling time ttun to be distinguishable is 




(4. 17) 



1 + lOfcol 

koa 



(4. 18) 



ttun » ^/M\ako\/ko. 



(4. 19) 
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We note that for a parity symmetric potential akg = —Mdhg/kQ, hence the 
condition becomes ttun » M dko/ko- For the long square barrier, this implies 
that 

-Y? d 

« 1. (4. 20) 

This condition can only be satisfied if jko/^o << 1. This is inadmissible, because 
the expansion (|3. 17p is not adequate in this regime. Hence, for the long square 
barrier the operational definition of the tunneling time is not meaningful. On 
the other hand, there is no problem in the short barrier limit {d ^ 0). 



4.3.3 The dependence on L 

Finally, we comment on the assumption that L » d. The consideration of a 
detector at a macroscopic distance away from the barrier region greatly simplifies 
our results: it leads to an expression for the tunneling time, which essentially 
coincides with the results of the asymptotic analysis of the wave packets. This 
assumption enters at two steps. First, in the construction of the POVM, we 
assume that L is sufRciently removed from the barrier region, so that the value 
for the particle's propagator sd x — L can be substituted by the corresponding 
value for the free particle. This condition is satisfied exactly if the corresponding 
Hamiltonian has no (generalized) eigenstates with an asymptotic behavior that 
does not correspond to that of a free particle (e.g. negative energy states). 
This is the case we considered in this section. Hence, the only place where 
the assumption of large L enters in a non-trivial way in the construction, is 
when we smear the probability distribution in order to remove the contribution 
of the terms oscillating as e**^"^. This implies that (at least formally), the 
expression (|4. 4p for the parameter Xkg is valid for all values of L such that 
the first condition stated above holds. We therefore obtain an expression for 
the tunneling time, even if the detector is located near the tunneling region. 
Clearly, this will have a very sensitive dependence on L, because the presence 
of the detector close to the barrier affects the configuration of the system. Note 
however that this result is rather formal, since it involves the idealization of the 
detection process by the crossing of the sharply defined surface x = L. In a 
realistic treatment the detailed physics of the detector are expected to influence 
the tunnelling time. 

For example, for a parity symmetric potential {Sk = 0), we obtain the 
following expression for the parameter Xkg , 

. rkoi2L + 9l)[l + cos(2fcoL + 9^,)] + r',^ sin(2fcoL + Ok,) 

where we wrote = r^.e*^'' and the prime denotes differentiation with respect 
to k. 
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5 A POVM for the tunneling time through se- 
quential measurements 

We saw in the last section that the determination of tunnehng time through the 
time-of-arrival probabihty is only meaningful for a specific class of initial states, 
because the delay time is not a proper random variable on the sample space 
of the POVM. It depends on the particle's initial momentum (and position) 
and as such it cannot be inferred unless both the initial state and the detection 
probability have very sharp maxima. 

However, this problem can be alleviated if we make a change in the ex- 
perimental set-up, namely if we consider that a measurement of momentum 
takes place before any recording of the time-of-arrival. In effect, if one consid- 
ers sequential measurements, it is possible to construct a POVM for which the 
tunneling time is a genuine random variable and no mixed classical-quantum 
arguments are needed for its identification. 

The procedure is the following. Let Q{x, k) be a POVM for unsharp phase 
space measurements. Let us also assume that the corresponding device is placed 
at the left-hand-side of the barrier; we perform an unsharp phase space measure- 
ment to any particle that moves towards the barrier that allows us to determine 
unsharp values for its position q and momentum p. The measurement is as- 
sumed to be non-destructive, hence the particles continue their motion, some of 
them cross the barrier and they are detected at distance L away. In other words, 
we have a sequential measurement: first an unsharp phase space measurement 
and then a time-of-arrival measurement. For each particle, the outcomes of this 
sequential measurement is encoded in the three numbers (x, k, t) that span a 
sample space f2. 

The key point is that from the knowledge of Q and II (the time-of-arrival 
POVM), it is possible to construct a POVM E on fl. The procedure is standard, 
see [24] for a detailed analysis. The POVM E consists of the positive operators 



E{t,x,k) = ^JQ{x,k)Il{t)^JQ{x,k), (5. 1) 

and of the positive operator 

E{N,x,k) = ^/Q{x,k)fl{N)/0{x,k), (5.2) 

that corresponds to a phase space measurement and then no detection. By 
construction it satisfies 

Jdxdk2TT(^J^ dtE{t,x,k) + E{N,x,k)^ =1. (5.3) 

For an initial state po, the joint probability density on the sample space fl 
is given by 

P{t, X, k) = Tr {poE{t, X, k)^ . (5. 4) 
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The key benefit in the consideration of such a POVM is that the delay-time 



and the tunneling time 



M(L -x) , , 

td^t ^ (5.5) 



ttun =td-\ ; — , (5. 6) 



are both random variables on the sample space Q. Hence, it is possible to 
define a POVM on the space in which they take values. We will do so after we 
construct expHcitly the POVM E. 

We consider POVMs for the unsharp phase-space measurements of the form 

— ^;^f{x - xo,k~ ko)\xo,ko){xo,ko\, (5. 7) 

where |a;o, kg) is the coherent state (|3. 12p . and / is a positive-valued function 
that determines the phase space resolution of the apparatus. Since f ^^^^Il{x,k) - 
1, it is necessary that the function / satisfies 

f(x,k) = l. (5.8) 

zvr 

The minimum resolution measurements correspond to f{x, k) = 27r6{x)d{k), 
in which case Q{x,k) = \xk){xk\. For simplicity, we will consider minimum 
resolution measurements in what follows. 

We obtain the following probability density on f2 

P{t,x,k) ^ {xk\p\xk){xk\fl{t)\xk). (5. 9) 

We note that {xk\Il{t)\xk) equals the probability density p{t) of Eq. (|3. 22p . 
We write this as Px,k{t), in order to express its dependence on the initial state 
[k = ko and x = in Eq. p. 22p ]. We then obtain 

P{t,x,k) = {xk\p\xk)px,k{t). (5. 10) 

We next change variables in (|5. 10[) from t to the delay time td- We note that 
on the full sample space, the relation between td and t is not one-to-one. First, 
the random variable td takes values in the whole real axis, while t only on the 
positive real axis. It is therefore convenient to define the probability P{t, x, k) 
for t running to all reals. This involves defining Px.k{t) for all t S R; we saw in 
Sec. 2 that this is obtained by doubhng the values of Px,k{t) for t e [0, oo). With 
t defined over all reals, we note that for each value of t, one obtains the same 
value for td twice, since is the same at points (t, x,p) and {t, 2L — x, —p). We 
perform the change of variables taking the facts above into account, and then 
we integrate over x and k, in order to obtain a marginal probability distribution 
over td 

Pd{td) = 4 / '^{xk\p\xk)pxM{td + ^^)- (5- 11) 
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The same procedure leads to a marginal probability distribution for the tunnel- 
ing time 

T-, / \ A f dxdk , , I ^1 , , L — X + dk . 

^tuniUun) ^ 4: / {xk\p\xk)Px^k{ttun ^ ^ )• (5. 12) 

The two equations above are completely general, and they hold without any 
approximations. They simplify significantly if we assume that for all values of k 
in the support of the initial state, the following two conditions hold: (i) Px,k{t) 
is appreciably different from zero only for times t such that t^a"^ /M << 1, and 
(ii) (T^fc << 1. The dependence on x of Px^k is then absorbed in the definition 
of the variable td^ and we obtain 



k\ 



Pd{td) = / dk{k\po\k)\Bk\^Ll exp <^ - 



M 



2fcV^ 



td- 



(5. 13) 



Similarly, 



nu„(W) = VW / dfc(fc|po|fc)|Bfc|2U exp. 



2fc2a_2 



M(Afe + dk) 



>5.. 



Note that neither Pd nor Ptun are normalized to unity. The delay and tunneling 
times are only defined for the fraction of the ensemble that corresponds to 
particles that have crossed the barrier. To normalize it, we have to divide by 
the probability corresponding to the detected particles 1 — Tr {j>qE{N)^ . 

Hence, we have constructed a positive definite probability density for the 
delay and the tunneling times, which is valid for an arbitrary initial state (with 
the restriction that its position support lies on the left side of the barrier) . This 
probability is definable in the context of a sequential measurement: there is no 
other way to define these times as quantum observables otherwise: the definition 
in Sec. 4 involved a mixture of quantum mechanics and classical argumentation 
and was only meaningful for a specific class of initial states. We have to keep 
in mind though that the experimental set-up for which these probabilities are 
valid involves keeping track of the phase space properties of individual particles 
and then comparing them with the registered arrival time. It requires relatively 
precise measurements at a microscopic scale, and it cannot be implemented 
when working with particle beams. 

We should also note that both probabilities Pd and Ptun are contextual, i.e. 
they depend strongly on specific features of the apparatus that performs the 
phase space sampling. They both have a strong dependence on the parameter 
cr, which defines the family of coherent states: in the present context cr is the 
inherent uncertainty in the specification of momenturr|j. At the limit ct ^ 0, 
both (|5. 13p and (|5. 14p vanish. There is, however, a limit in which the results 
become cr- independent. If the initial state has support on values of fc, such that 



^For the contextuality of sequential measurements, see the extended discussion in 
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the mean of the Gaussian in either probabihty density is much larger than its 
spread, then we can approximate it by a delta function. This condition implies 

<7\\k\ » 1, (5. 15) 

for ([5~T^ and 

a{Xk+dk) » 1 (5. 16) 

for (I5~T4)1 . 

At these regimes, we obtain 

Pd{td) ^27tJ dk{k\po\k)\Bk\H{td ~ (5. 17) 

PtuniUun) ^21, j dk{k\p^\k)\Bk\H{U^^ - (5. 18) 

In other words, the values of Pd{td) and of Ptun{itun) are determined by the 
value of the probability distribution of the initial's state momentum at values of 
k that are solutions of the algebraic equations td = ^^^^ and ttun — ^'^('^^+'^''^ 
respectively. These expressions for the probability distribution are independent 
of the detailed characteristics of the phase space POVM: they only depend on 
the initial state and on the characteristics of the potential^. They can therefore 
be considered as ideal distributions of delay and tunneling times respectively 
that exhibit little sensitivity to the measurement scheme employed for their 
determination. 

We can further simplify the expressions for Pd and Ptun using the estimation 
dOl) for Bko-. 

Pd{td) = j dk{k\Mk)\Tk\H{td - ^), (5. 19) 

Ptun{ttun) = j dk{k\po\k)\n\H{Uun " M^^lM). (5. 20) 

In effect, the probability for td and ttun are defined from the corresponding 
values of the momentum distribution weighted by the transmission probability. 
Defining the functions Fd{k) :— and Ftun{k) '■— :^ll^i±f^i)^ -^e see that the 
probabilities ()5. 19115. 20|) are obtainable from the operators Td = Fd{p) and 
Ttun = Ftunip) (p is the momentum operator) when these act on the state 

Pcross = J dk\Tk\''PkPoPk, Pk = \k){k\ (5.21) 

that describes the sub-ensemble of particles that have crossed the barrier. One 
could therefore call Td and Ttun time-delay and tunneling-time operators re- 
spectivelj0. 



^Recall that by virtue of smearing the value of L, there is no L-dependence in Bf^ ; hence 
the marginal probability distributions are also L-independent. 

■^There is an ambiguity in their definition at fc = 0. However, this does not affect the 
probabilities l|5. 19| l and l|5. 20|l . because ITfe^ol = 0. 
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We end this section, by examining the domain of vahdity of conditions (|5. 16p 
and (|5. 15[) for the square potential barrier. At the large barrier limit, they read 



a\-dA 1 >> 1 

7fe 

cr/7fc » 1. 



(5. 22) 



(5. 23) 



They are satisfied if the position a^^ spread of the coherent states is much 
smaller than the effective lengths corresponding to delay and tunneling time 
respectively. For the delta function barrier, these conditions imply 



which is only possible if n is extremely small (a rather unphysical case). 

We see therefore that the ideal probability distributions (|5. 17p and (|5. ISp 
can only be obtained if the initial phase space measurement has a resolution 
for position substantially smaller than the dimensions of the barrier. This is 
a type of measurement that is not explicitly forbidden by quantum mechanics, 
but clearly it would be extremely difficult to achieve in practice. 

6 Conclusions 

We reformulated tunneling as a problem in the determination of probability for 
the time-of-arrival. This allowed us to identify the classic Bohm-Wigner time as 
the most suitable measure for the tunneling time. However, this identification 
only holds for a specific class of initial states and potentials; in other regimes, 
there is no operational definition of the concept. There is one way to go around 
this problem by considering a sequential measurement set-up: we first measure 
the phase space properties of the particles (before they attempt to cross the 
barrier) and then we determine their times-of-arrival. In this context, it is 
possible to construct a probability measure for the tunneling time that is valid 
for all initial states. 

The key feature of our construction is that there is neither interpretational 
nor probabilistic ambiguity. The probabilities we derive are obtained through 
a POVM, hence they are always positive and they respect the convexity of 
the space of quantum states. The interpretation of these objects is concretely 
operational, in the sense that it is tied to the statistics for the measurement of 
particles' arrival times. Tunneling time is solely defined in terms of the statistics 
of measurement outcomes. 

In another paper [2], the POVM wc constructed here will be employed for 
the study of the decay probability of unstable states through tunneling. 
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